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We study the thermoelectric properties of three-dimensional topological insulators with many 
holes (or pores) in the bulk. We show that at high density of these holes the thermoelectric figure of 
merit ZT can be large due to the contribution of the conducting surfaces and the suppressed phonon 
thermal conductivity. The maximum efficiency can be tuned by an induced gap in the surface states 
dispersion through tunneling or external magnetic fields. The large values of ZT, much higher 
than unity for reasonable parameters, make this system a strong candidate for applications in heat 
management of nanodevices, especially at low temperatures. 



Efficient heat management in nanodevices and energy 
conversion from wasted heat have been some of the key 
driving motivations for the search of new thermoelectric 
materials pl-[T4l. Their efficiency rests on the ability of 
a material to have a low thermal conductivity and high 
thermopower and electric conductivity. The large ther- 
mopower requires steep dependence of the electronic den- 
sity of states on energy, which can be achieved by having 
the chemical potential close to the bottom of a band. 
The relatively high conductivity demands the gap to be 
low. These combined requirements point to semiconduc- 
tors with heavy elements as the best candidates for the 
thermoelectric materials. 

A typical band gap for an efficient thermoelectric de- 
vice at room temperature should be of the order of few 
hundred meV. Such a low gap can appear in the mate- 
rials either accidentally or due to a large enough spin- 
orbit coupling that leads to the inversion of the band 
structure, thus providing a band-gap of the order of the 
spin-orbit coupling energy. This latter class of materi- 
als, besides being highlighted by their high thermoelec- 
tric efficiency, has been identified as having an insulat- 
ing bulk and topologically protected conducting surfaces 
with Dirac-like band structure p!5HT8] : hence their name, 
topological insulators (TI). 

The bulk thermoelectric properties of this family of 
semiconductors have been extensively studied both the- 
oretically [6'-'8^ and experimentally [TQHT3]. The rather 
high thermoelectric efficiency found in these materials is 
the property of the band structure and the heavy atomic 
masses, and is unrelated to the contribution of the pro- 
tected surface states — the main focus of this paper. 
As with most thermoelectric materials, the main stum- 
bling block for increasing the thermoelectric efficiency is 
the high phonon thermal conductivity. There have been 
many attempts to decrease it while leaving the electric 
conductivity intact [19, 20 , e.g., by introducing disor- 
der in the material that suppresses phonon transport 
(phononic glass). The bulk electronic states are how- 
ever often affected by this disorder as well and it is very 
difficult to reduce thermal conductivity without also re- 
ducing thermopower and electric conductivity [Ij. 

Here we propose a new approach to increase ther- 
moelectric efficiency by combining the novel properties 
of the TI materials, i.e., their protected surface states. 



FIG. 1. A sample made of a holey 3D topological insulator. 
These holes increase the surface to bulk ratio and reduce its 
phonon thermal conductivity. The inset shows two such holes 
with the Dirac like metallic surface states. 



with methods to reduce thermal conductivity that can 
enhance the contribution to the conductivity and ther- 
mopower from these topological surface states by increas- 
ing the surface to volume ratio. 

The thermopower efficiency is characterized by the 
dimensionless thermoelectric figure of merit ZT = 
(jS'^T/ n^ where a, S^ n^ and T are the electric conduc- 
tivity, thermopower, thermal conductivity, and tempera- 
ture, respectively. This quantity corresponds to the ratio 
of the output power to the rate of the heat power con- 
sumption under an applied temperature gradient. 

The high surface to bulk ratio can be achieved by 
nanostructuring 3D topological insulators with many 
holes in the transport direction. The surface of each hole 
behaves as a 2D conductor. The effect of the holes is 
threefold: 1) to suppress phonon thermal conductivity - 
the strongest limiting factor for ZT, 2) to increase con- 
ductivity due to a high contribution of conducting sur- 
face, and 3) to allow for the surface states to tunnel and 
create a smaller subgap to increase ZT further. A sketch 
of the sample made of a 3D topological insulator with a 
high density of holes is depicted in Fig. [l] Even at high 
density of these holes, the topologically protected surface 
states remain robust. It has been shown that these states 
survive in the films as thin as of three or more quintuple 



layers [2T] . 

Reducing thermal conductivity caused by the increased 
phonon scattering off grain boundaries and defects has 
been experimentally confirmed recently. For example, 
in BiSbTe bulk alloys higher ZT has been observed in 
experiments with nanocrystalline bulk samples made by 
pressing nanopowders at high temperatures [19 . Fur- 
thermore, high density of periodic nanoscale holes in sil- 
icon thin films has been shown to reduce the thermal 
conductivity by two orders of magnitude [20^ . 

The unique feature of the topological insulators is that, 
although they have a bulk energy gap of 2Ao, they 
also have ungaped surface states. These propagating 
states are confined to the close proximity of the sur- 
faces. Their existence is protected by the topology of 
the bulk band structure. These states have cone-like 2D 
Dirac spectrum, E = ^vh\k\^ where v is the constant 
Dirac electron velocity [22 . Application of a magnetic 
field or hybridization of these states due to close prox- 
imity of two surfaces can induced a controllable Dirac 
subgap A [23^, ^24^. Then the spectrum takes the form 
E = ^\/v'^h?k?' + A^. Generally this subgap is much 
smaller than the bulk gap, A <C Aq. For 3D TIs 2Ao is 
roughly from 0.15 eV (for Bis Tea) to 0.3 eV (for BisSea). 

To study the transport properties of these materi- 
als we use linear response theory [25l [26]. The elec- 
tric (j^) and thermal {j^) currents are given by lin- 
ear combinations of the chemical potential and tem- 
perature gradients: j^/e = LqV ji + I/i(VT)/T and 
j^ = — LiV/i — L2(VT)/T, where e is the electron charge. 
Using Onsager relations, one can find from these equa- 
tions the electrical conductivity a = e^Lo, thermopower 
S = — Li/(eTLo), and electronic thermal conductivity 
Ke = {L0L2 - Ll)/{TLo). The figure of merit, ZT, can 
then be represented in terms of these linear coefficients 
as 
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Here K:ph is the phonon contribution to the thermal con- 
ductivity in the bulk (phonon contribution to the ther- 
mal conductivity for the surface TI states is much smaller 
than that in the bulk). In Eq. (fTj) it is assumed that the 
transport coefficients have bulk and surface contributions 

Ln = Lb,n^Ls,n/V, whcrC V = (A-En ^^n)/ En ^^^n 

is the factor related to surface/bulk ratio (porosity) and 
has dimension of length. The parameter V characterizes 
the average distance between the pores (holes) of the av- 
erage radius R, see Fig. |2](b). 

It follows from Eq. (IT]) that at small P, the contribution 
to ZT mostly comes from 2D surface states and in this 
limit ZT is given by 
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ZT2D is shown in the inset of Fig. [3] for A/(/cbT) = 3 
as a function of the chemical potential. The color plot of 
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FIG. 2. (a) A sample of 3D topological insulator with ran- 
domly distributed pores, (b) Top view of small part of the 
sample with the holes. V represents the average distance be- 
tween randomly distributed pores (or holes). 
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FIG. 3. The color plot of ZT only for the surface states 
as a function of the induced subgap A/{kBT) and chemical 
potential fi/{kBT) measured in units of temperature. The 
inset shows ZT calculated only for the surface states with the 
subgap A/{kBT) = 3. The solid golden line represents the 
maximum ZT at a fixed A. The black line corresponds to 



ZT2D as a function of both the induced subgap A/iksT) 
and chemical potential /i/{kBT) is shown in Fig. [3J For 
a fixed induced gap A/(/cbT) we plot the location of the 
maximum ZT achievable by tuning the chemical poten- 
tial (golden line in Fig. [3| . The very high values of ZT 
in reality are not reachable since any small contribution 
from the phonon thermal conductivity will reduce ZT. 

The phononic contribution to the thermal conductiv- 
ity for the holey sample can be characterized by the di- 
mensionless parameter Kph = 2hZphVh'^ / {rk^T'^) . Here 
the phonon thermal conductivity can be estimated to be 
tvph ~ lWm~ K~^ (as for Bi2Te3) and the average dis- 
tance between the holes reaching P '^ 10 nm. Taking the 
relaxation time r ~ 10~^^ s, at room temperatures we 
estimate Kp^ ~ 1 for these rather conservative values of 
parameters. In Fig. [4] we show ZT as a function of the in- 
duced subgap A/{kBT) and chemical potential ii/{kBT) 
for the phononic bulk contribution characterized by the 
dimensionless parameter Kp}^ = 3. Although reduced 
considerably from its pure 2D value, ZT remains sub- 
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FIG. 4. ZT as a function of the induced subgap /\/{kBT) 
and chemical potential ii/{kBT) with dimensionless param- 
eter Kph — 3 characterizing the phonon contribution to the 
thermal conductivity from the bulk states. The solid orange 
line represents the maximum ZT at a fixed A. 



stantially larger than any value so far achieved in these 
materials and can be tuned significantly by changing the 
chemical potential, e. g., by gating the system. 

We also investigate the relation between the induced 
subgap A/{kBT) and the chemical potential lJimax/{^BT) 
to obtain the maximum value of ZT for different 
phononic bulk contributions. This result is shown in 
Fig. IsJ For small phonon thermal conductivity (small 
Kph) the ZT is maximal when the chemical potential is 
placed within the subgap, /i < A. In the limit of large 
Kph the maximum ZT is reached when the chemical po- 
tential and the gap match. This can be seen in Fig. |5J 
where the line corresponding to K^^ = 90 almost coin- 
cides with the dashed line A = jirnax for rather large 

l-J^max/ikBT). 

Finally, we point out that the large values of ZT ob- 
tained in the absence of phonon thermal conductivity 
shown in Fig. [3] is not the entire story. As at large 
A/{kBT) and ii/{kBT) these thermoelectrics are indeed 
very efficient but not effective - the power factor aS'^^ 
describing how much power one can produce, is small 
in this limit. This is shown in Fig. [6] as a plot of the 
maximal ZT as a function of the power factor. The di- 
mensionless power output L\/Lq cx cfS'^ as a function 
of the induced subgap A/{kBT) and chemical potential 
li/{kBT) is shown in the inset of Fig. [6J 

Within these materials, there are extensive competing 
mechanisms that contribute to the thermal transport and 
thermoelectric efficiency. In order to distinguish that a 
particular increase of ZT is associated with the topolog- 
ical protected states, we propose a simple measurement 
of the transport coefficients with and without a magnetic 
field applied, parallel to the transport direction, to in- 
crease the induced gap of the surface states in the energy 
spectrum. The effect of the magnetic field will be more 
pronounced when the temperature is smaller than the in- 
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FIG. 5. Induced subgap A/{kBT) as a function of chemical 
potential iimax/{kBT) for the maximum value of ZT. The 
dimensionless parameter Kph characterizing the contribution 
of the phonon thermal conductivity is taken to be 0; 1; 3; 9; 
and 90. The dashed line iimax = A is shown for comparison. 
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FIG. 6. Maximum ZT as a function of the dimensionless 
power factor, L\/Lq oc aS^ . The inset shows the power fac- 
tor as a function of the induced gap A/(kBT) and chemical 
potential ii/(kBT) measured in units of temperature. 



duced gap. For fields of ~ 10 Tesla a gap of the order 
of several meV is expected in materials such as Bi2Se3 

[271 Eg. 

It is important to note that our consideration is not 
limited to periodic holey structures. It is also applicable 
to porous materials with random size and location of the 
holes (pores) as shown in Fig. [2|a). In this case some of 
the surface states can be situated too close to each other 
so that the back-scattering in them will be increased, 
meanwhile the phonon part of thermal conductivity will 
be reduced due to stronger trapping of phonons in chaotic 
structures. Therefore, ZT in porous materials can be 
comparable to those with the periodically placed holes. 



METHODS 

To estimate the surface contribution of TI to the trans- 
port coefficients i„ we assume the bands to be Dirac-hke 
with a subgap A and use Boltzmann equation in the re- 
laxation time approximation, 
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Here the sum is over the upper and lower bands, i = ±1, 
and f{E) = df/dE with / = l/(e(^-^)/(^^^) + 1) being 
the Fermi distribution function. Then taking relaxation 
time r to be independent of energy, we find 



potential is far below the bottom of the conduction band, 
(Ao — /j.)/{kBT) :^ 1, the contribution from the bulk 
to the transport coefficients is exponentially suppressed, 
and can be safely neglected [29]. 
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The same applies for the valence band contribution when 
the chemical potential is far away from the band edge. 

Then the thermoelectric transport is dominated by the 
surface states and the only sensible contribution from the 
bulk is to the phonon thermal conductivity hZph^ so that 
the figure of merit takes the form 
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where h is Planck constant, A = A/{kBT), and /i 

To estimate the bulk contribution to Ln we assume 
the bands to be parabolic and find in the relaxation time 
approximation for the conduction band, 

L,,„ = -rj^ D{E)^-^j f{E){E-^^rdE, (5) 

where D{E) is the density of states. When the chemical 



This expression was used to obtain the results shown in 
Figs.[4]and|5j 
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